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Cross-Reference to Related Applications 

This patent is a continuation-in-part of International Application 
PCT/CA02/0069 designating the United States. 

Technical Field 

This patent relates to signal processing techniques and systems and, in 
particular, to a method of extraction of sinusoidal signals of time-varying nature and an 
apparatus built on the basis of such a method. 

Background 

The effect of periodic phenomena, registered as signals, is usually studied 
using common signal analysis tools such as Fourier analysis. A given signal, as long as it 
manifests some periodicity, can be thought of as being composed of a series of single 
sinusoidal components or, more commonly named, tones. Frequency domain characterization 
of signals amounts to the identification of individual constituting components of a given 
signal and their corresponding values of amplitude and phase. In the case of nonstationary 
signals, all characteristics of the constituting components, including frequency, may vary 
over time. The main shortcoming of Fourier-based methods is their inherent fixed-frequency 
assumption which limits their applicability to the majority of real signals, the frequency 
characteristics of which may vary with time. There have been numerous attempts to devise 
adaptation mechanisms to be incorporated into signal analysis tools to render them useful for 
analysis of quasi-periodic signals. Linear adaptive filtering is an example of such methods 
reported so far with its success and shortcomings. 

Aside from the time-varying quality of real signals, very often signals are 
buried under noise and disturbances and may thus be severely distorted. Often, a useful signal 



analysis tool loses its efficiency when it is applied to signals affected by noise and 
disturbances. Thus, it is often necessary to recover the signal itself out of the background 
noise, especially under nonstationary conditions. Extraction of the signal itself, and not just 
its characteristics, is of particular interest in applications where synchronization matters; i.e. 
5 where the total phase information of a signal is important. In such cases, a single sinusoidal 
component of a given signal, or the totality of a number of such sinusoids, is to be extracted, 
or equivalently a desirable noise-free synchronized signal is to be synthesized. Synthesization 
of signals synchronous with a jgiven reference signal finds applications in more areas than 
those dealing with extraction of signals out of noise and may very well include those 

1 0 applications in which phase-locked loop (PLL) circuits and systems are employed. k 
..... Considering the inadequacy, of the performance of the available standard, 
signal analysis tools, such as Fourier-based techniques, adaptive filters and PLLs, in 
extracting sinusoids of time- varying nature buried under noise in a unified way, it is not 
difficult to explain the existing diversity of the methods, each designed to tackle a specific 

15 type of problem. 

. Signal processing techniques and systems, and in particular methods of 
extraction of a desired signal within a noisy environment find significant application in v 
audiometry, i.e., the measuring of how well someone can hear. Thus improvements in 
techniques, systems and methods for the extraction of a signal from a noisy environment can 

20 lead to improvements in audiometry. 

Conventional audiometry is performed by having a subject respond to acoustic 
stimuli by pressing a button, saying "yes", or repeating words that may be presented in the 
stimulus. These tests are subjective in nature. Audiometry allows an audiologist to 
determine the auditory threshold of the subject, which is defined as the lowest intensity at 

25 which a sound can be heard. The audiologist evaluates the auditory threshold of a subject by 
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using a stimulus that most commonly consists of a pure tone. The stimulus is presented via 
earphones, headphones, free field speakers or bone conduction transducers. The results are 
presented as an audiogram which shows auditory thresholds for tones of different 
frequencies. The audiogram is helpful for diagnosing the type of hearing loss a subject may 
5 have. The audiogram can also be used to fit a hearing aid and adjust the level of 
amplification of the hearing aid for subjects who require hearing aids. 

Conventional audiometry cannot be performed if the subject is an infant, 
young child or cognitively impaired adult. In these cases, obj ective tests of hearing are 
necessary in which the subject does not have to make a conscious response. Objective 

10 audiometry is essential for detecting hearing impairment in infants or elderly patients as; well 
as for evaluating functional hearing losses. Furthermore^ few objective tests have been 
developed for supra-threshold tests of speech, frequency, or intensity discrimination. 

One form of objective audiometry uses. auditory evoked potentials. Auditory: 
evoked potential testing consists of presenting the subject with an acoustic stimulus and 

15 simultaneously or concurrently sensing (i.e. recording) potentials from the subject. The nv 
sensed potentials are the subject's , electroencephalogram (EEG) which contain the subject's , 
response to the stimulus if the subject's auditory system has processed the stimulus. These 
potentials are analyzed to determine whether they contain a response to the acoustic stimulus 
or not. Auditory evoked potentials have been used to determine auditory thresholds and 

20 hearing at specific frequencies. 

One particular class of auditory evoked potentials is steady-state auditory 

evoked potentials (SS AEPs). The stimulus for the SSAEP consists of a carrier signal, which 

- 

is usually a sinusoid, that is amplitude modulated by a modulation signal which is also 
usually a sinusoid. The SSAEP stimulus is presented to the subject while simultaneously 
25 recording the subject's EEG. If the auditory system of the subject responded to the SSAEP 
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stimulus, then a corresponding steady-state sinusoidal signal should exist in the recorded 
EEC The signal should have a frequency that is the same as the frequency of the modulation 
signal (i.e. modulation frequency). The presence of such a corresponding signal in the EEG 
is indicative of a response to the SSAEP stimulus. Alternatively, the phase of the carrier 
5 signal may be frequency modulated instead of or in addition to amplitude modulation to 
create the SSAEP stimulus. A system and method of SSAEP audiometry is described in 
United States Patent No. 6,602,202, the disclosure of which is hereby expressly incorporated 
by reference. - 

However, objective audiometry employing SSAEP testing is time-consuming 

10 because the amplitude of the SSAEP response is quite small compared to the background v 
noise^ which is.the subject's ongoing brain activity (i.e. EEG) while the test is being 
conducted. The SSAEP response thus has a small signal-to-noise ratio (SNR) which makes it 
difficult to detect the SSAEP response in a short time period. Thus, improvements in v 
techniques, systems and methods of detecting the SSAEP response may lead to 

15 improvements in SSAEP audiometry. 

Another form of objective audiometry is otoacoustic emission (OAE) 
audiometry. One form of OAE is distortion product otoacoustic emission adiometry 
(DPOAE). Distortion product otoacoustic emissions (DPOAEs) are very low level stimulated 
acoustic responses to two pure tones presented to the ear canal. DPOAE measurement 

20 provides an objective non-invasive measure of peripheral auditory function and is used for 
hearing assessment. DPOAE screening is becoming a standard.clinical practice to predict 
potential sensorineural hearing loss especially in newborns. 

DPOAEs have been recognized for a number years. But DPOAE measurement 
remains difficult because of the challenging nature of the signal processing task. In this type 

25 of otoacoustic test, two pure tones with frequencies /, and f 2 are presented to the cochlea. 



For best results, f 2 is usually chosen as 1.2 /, . Due to the non-linearity of the ear, a very low 
level of distortion product of frequency 2 f l -f 2 is generated in normal ears. The level of such 
DPOAE signal is a measure of functionality of the ear. Estimation of such a weak signal 
buried under two strong artifacts in a potentially noisy background is a challenging signal 
processing problem. 

Conventionally, fast Fourier transform (FFT) is used as the main signal 
processing tool to estimate the level of DPOAE signals. Application of FFT in this problem 
has a number of shortcomings among .which long measurement time is most pronounced. 
Such long measurement time is usually required for acquisition of more data which, when 
averaged, reduce the overall background noise effect. Unreliability of the measurements is 
another problem of FFT-based methods and is a direct result of the sensitivity of the FFT- 
based methods to the background noise. In addition to the need to increase the measurement 
time, the tests are usually required to be conducted in low noise environments such as sound- 
proof booths. 

In an attempt to devise high performance DPOAE estimation techniques, 
linear adaptive signal processing techniques have been employed. Such techniques generally 
offer better performance in terms of measurement time, which may be interpreted as higher 
noise immunity of adaptive techniques compared to FFT; however, the need for sound-proof 
examination rooms is not obviated with such techniques. 

Brief Description of the Drawings 

In the accompanying drawings like reference numbers denote like 
components, brief description of each of which is herewith given. 

Figure 1 shows a general structural block diagram of a system embodiment 
described herein; 
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Figure 2 illustrates, by way of example, the performance of a system or 
method according to the described embodiments described herein in extracting a sinusoidal 
signal; 

Figure 3 illustrates, by way of example, the performance of a system or 
5 method according to the described embodiments described herein in tracking a step change in 
the amplitude of the input signal; 

Figure 4 illustrates, by way of example, the performance of a system or 
method according to the described embodiments described herein in tracking a step change in 
the frequency of the input signal; * 
10 v ; Figure 5 illustrates, by way of example, the, effect of the presence of noise, 

: Figure 6 depicts a functional block diagram of a system according to another 

of the described embodiments; 

Figure 7 depicts a functional block diagram of a system according to another; 
of the described embodiments adapted as a phase-locked loop (PLL); 
15 Figure 8 depicts a functional block diagram of a system according to another, 

of the described embodiments incorporating a filter to enhance its performance; 

Figure 9a depicts a functional block diagram of a system according to another 
of the described embodiments incorporating a plurality of core elements in parallel 
combination; 

20 Figure 9b depicts a functional block diagram of a system according to another 

of the described embodiments incorporating a plurality of core elements in cascade 
combination 

Figure 10 depicts a functional block diagram of a system according to another 
of the described embodiments adapted for the elimination of undesirable components; 
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Figure 1 1 depicts a function block diagram of a audiometer adaptable to 
operate in accordance with one or more of the described embodiments; 

Figure 12 is a detailed block diagram of the probe component and the data 
acquisition component of the audiometer depicted in Figure 1 1 ; 
5 _ Figure 13 is a functional block diagram illustrating a software architecture that 

may be used in the audiometer depicted in Figure 11; 
' s -\. Figure 14 is a functional block diagram of a signal processing scheme that 

may be used:in the audiometer depicted in Figure 1 1 for providing DPOAE audiometery; 

Figure 15 illustrates, by way of example, performance of an audiometer 
TO configured according to one of the described embodiments for the extraction of DPOAE; 

. Figure 16 illustrates, by way of example, performance of an audiometer 
configured according to one of the described embodiments for the extraction of DPOAE > 
responsive to a step change; 

Figure 17 illustrates, by way of example, performance of an audiometer 
1 5 configured according to one of the described embodiments for the extraction of DPOAE fori 

multiple frequencies; 

■ < • 

: Figures 18a and 18b illustrate, by way of example, performance of an 
audiometer configured according to one of the described embodiments for the extraction of 
DPOAE in a noisy environment, wherein Figure 18b is an expanded portion of Figure 18a; 
20 and 

Figures 19a and 19b illustrate, by way of example, performance of an 
audiometer configured according to one of the described embodiments for the extraction of 
DPOAE in a noisy environment, wherein Figure 19b is an expanded portion of Figure 19a. 
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Detailed Description 

This patent describes embodiments of apparatus, systems, algorithms and 
methods of extraction and estimation of the parameters of individual constituting sinusoids 
present in a given signal. The described embodiments offer an alternative to standard Fourier 
5 transform or linear adaptive filtering as regards signal analysis, and to phase-locked loops as 
regards signal synthesis while maintaining a structural simplicity and high speed of :* H 
convergence both comparable with those of Fourier-based techniques. Unlike Fourier-based \ 
techniques, once the individual constituting components of a given signal are extracted, time 
variations in the characteristics of the components -are registered and tracked. ; 

10 : An apparatus, system, algorithm or method in accordance with one or more of 

the- described embodiments may provide for the extraction of a more or less specified desiredv 
sinusoidal component of a given signal, characteristics of which such as its amplitude, phase 
and frequency may vary with time. Accordingly, synthesization of a signal synchronous with 
a desired component of a given reference time-varying signal is possible, thereby rendering a 

15 signal synthesis tool. .. . 

An apparatus, system, algorithm or method in accordance with one or more of ; 
the described embodiments may provide for the estimation of time-varying parameters of a 
more or less specified sinusoidal component of a given signal such as its amplitude, constant 
and total phase, and frequency. Accordingly, detection and measurement of characteristics of 

20 constituting components of a given signal is possible, thereby rendering it a time-frequency 
analysis tool. 

An apparatus, system, algorithm or method in accordance with one or more of 
the described embodiments may provide for the extraction of a more or less specified 
sinusoidal component of a given signal which may be highly polluted by noise and may be 
25 distorted by various external disturbances. Advantageously, a noise- free signal synchronous 
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with a desired component of a given time-varying signal which may be polluted by noise may 
be detected, thus featuring it as a noise elimination/reduction technique. Also 
advantageously, a high degree of noise immunity is possible, thus a tool suitable for analysis 
of intensely noisy time-varying signals is possible. 
5 An apparatus, system, algorithm or method in accordance with one or more of 

the described embodiments may provide structural simplicity so that the computational needs 
are comparable with those of simple Fourier-based techniques. 

"v Ail apparatus, system, algorithm or method in accordance with one or more of 
the described embodiments may provide for enhanced efficiency in terms of convergence 

10 time comparable; with that of /Fourier-based techniques. ^ r 

4 The apparatus^ systems, algorithms and methods according to the described 
embodiments relate to; signal processing techniques and systems and, in particular, to a j; 
method of extraction of sinusoidal signals of time-varying nature. An apparatus built on the 
basis of such a method is capable. of extracting a specified single sinusoidal component of an 

15 input signal, potentially containing other components and noise, and tracking variations of the 
i amplitude, phase and frequency of such a sinusoid over time. Various apparatus, systeriis^ 
algorithms and methods according to the described embodiments *ire generally referred to 
herein as a signal analysis/synthesis tool or simply "tool" or "core unit," which may be. used 
in isolation or as the fundamental building'block of single-core or multi-core systems. This 

20 tool or core unit therefore finds application in diverse areas of engineering arid science, 

ranging from signal detection, extraction, measurement and synthesis in engineering to time- 
frequency analysis of mechanically generated signals in science. Some examples of its 
applications are the detection of dual-tone multi-frequency (DTMF) signals in telephony, 
signal recovery in noisy biopotential signals, audiometry and extraction and measurement of 

25 power line signals. 
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The instant disclosure is provided to further explain in an enabling fashion the 
best modes of making and using various embodiments in accordance with the present 
invention. The disclosure is further offered to enhance an understanding and appreciation for 
the inventive principles and advantages thereof, rather than to limit in any manner the 
5 invention. The invention is defined solely by the appended claims including any amendments 
made during the pendency of this application and all equivalents of those claims, as issued. i. 7 ," 

• It is further understood that the use of relational terms, if any, such as first and 
second, top and bottom; and the like are used solely to distinguish one from another entity or 
action without necessarily requiring or implying any actual such relationship or order 

10 between such entities or actions. , : a i&; 

>;:; ; Much of the inventive functionality and many of the inventive principles are 
best implemented with or in software programs or instructions and integrated circuits (ICs);, > 
such as application specific ICs. . It is expected that one of ordinary skill, notwithstanding /> ; 
possibly significant effort and many design choices motivated by, for example, available 

15 time, current technology, and economic considerations, when guided by the concepts and 
principles disclosed herein will be readily capable of generating such software instructions^ 
and programs and ICs with minimal experimentation. Therefore, in the interest of brevity 
and minimization of any risk of obscuring the principles and concepts in accordance to the 
present invention, further discussion of such software and ICs, if any, will be limited to the 

20 essentials with respect to the principles and concepts of the preferred embodiments. 

*. According to the described embodiments, a tool for extraction of individual 
constituting sinusoidal components of a given signal may include a number of simple steps to 
provide estimated values of the amplitude, phase and frequency of a more or less specified 
component of a given signal as well as the synthesized desired component itself. Referring to 

25 Figure 1, let u(t) represent a signal which is registered by a sensor to indicate a natural 

-10- 



phenomenon or is artificially generated by the use of some equipment. It is usually taken to 
be a function of time t although it can equivalently be a function of any other variable, and 
may be of an electromagnetic nature, .such as voltage or current signals in electrical systems, 
or of a mechanical nature such as sound waves. Very often, signals exhibit some periodicity 
5 in which case they can be represented as a series of individual sinusoidal components. If a 
signal is not perfectly periodic and its characteristics happen to change over time, so will the 
parameters of constituting sinusoidal components. In such cases, it is desirable to extract 
individual. constituting sinusoidal components and track their variations over time. 

The input signal u(t) is used to obtain an output signal y(t) which is the desired 

10 constituting component of such a given signal. Such a synthesized output signal is subtracted 
from the given input signal by a subtraction operation 5. The outcome of such an operation is 
another signal e(t) which, by the very fact that it is the difference between the given input v 
signal and the synthesized desired component, is the totality of the undesirable components/ 
and the error incurred in the extraction process. This operation can be concisely formulated 

15 by the following equation: . > -y. 

:vr.;.: V e(t)=u(t) -y(t). < (1) 

;The instantaneous rate of change, or equivalently time derivative, of the 
amplitude of the desired component of the signal is taken to be proportional to the error 
signal e(t) and the sine of the total estimated phase of the desired component, and hence it is 

20 proportional to their product. Therefore, a multiplication operation 10, provides the product 
of the error signal e(t) and the sine of the total estimated phase of the desired component. 
Thus, a proportionality constant /u x 15, when multiplied by the outcome of the product 
operation 10, yields the instantaneous rate of change of the amplitude of the desired 
component which, when integrated by an integration operation 20, provides the estimated 

25 value of the amplitude of the desired component of the input signal. The value of the initial 
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condition of the integration 20 can be taken as zero in which case the tool is initialized from 
zero amplitude. It is easy to formulate the procedure just outlined for the estimation of 
instantaneous value of the amplitude of the desired component of a given input signal by the 
following equation: 

5 A = ju } e sin & (2) 

where the dot on top denotes the time derivative, and A and 0 denote the amplitude and 
phase, respectively- : 

Quite similarly, the instantaneous rate of change of the frequency of the desired 
component of the input signal may be taken to be proportional to the error signal e(t), the; 

10 . estimated value of amplitude A(t) and the cosine of the total estimated phase of the desired 

component, and hence it is proportional to their product. Therefore, a multiplication operation 
30 provides the product of the error signal e(t), the estimated value of the amplitude A(t) and 
the cosine of the total estimated phase of the desired component. Thus, a proportionality \ \ 
constant u 2 40, when multiplied by the outcome of the product operation 30, yields the r 

1 5 instantaneous rate of change of the frequency of the desired component which, when *; 
integrated by an integration operation 60, provides the estimated value of the frequency of the 
desired component of the- input signal. The value of the initial condition of the integration 60 
can either be taken as zero in which case the tool finds the sinusoidal component of the input 
signal whose frequency is closest to zero, or can be set to any desired initial frequency 
= a± . '. • 

Jo _ 

20 ln > in which case the tool finds the sinusoidal component of the input signal the 

frequency of which is closest to such a desired value^. This conveniently provides for 
specifying the desired component of the input signal the extraction of which is sought. The 
extracted component of the signal will therefore be the one angular frequency that is closest 
to the value of the initial condition co 0 50. It is easy to configure the tool for the estimation of 

25 the instantaneous value of the frequency of the desired component by the following equation: 

cb = ju 2 eAcos <p, (3) 
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where co represents the value of the angular frequency of the desired component of the input 

signal; frequency/in Hz is ^ Jr . 

The instantaneous rate of change of the total phase of the desired component 
of the input signal is the sum of the angular frequency co itself and a constant factor of its 
5 instantaneous rate of change. This can be formulated as: 

<j> = co + ju 3 a>, * (4) 

or with reference to the procedure of obtaining the estimation of frequency formulated by 
equation numbered (3), as: 

<fi = (D + /u 2 /j 3 eA cos <f>\ 

10 Thus, a combined constant factor ju 2 ju 3 35 provides the scaled product which 

is then added to the angular frequency by an addition operation 45 to yield the instantaneous 
rate of change of the total phase of the desired component of the signal which, when :j>h 
integrated by integration 55, yields the estimated value of the phase of the desired component 
of the input signal. The value of the initial condition of the integration 55 can be taken as zero 

15 in which case the tool is initialized from zero phase. The sine function 65 and the cosine 

function 70 generate the sine and cosine of the estimated total phase. The function of the sine 
operation 65 and cosine operation 70 can be interchanged without any effect oh the 
performance since this would only mean a different initial condition of integration 55. 

The desired component of the input signal is a single sinusoid having an 

20 amplitude of Aft) and a phase of ^ (t) and can therefore be generated by multiplying the 

amplitude and the sine of the phase of the desired component by a product operation 25. This 
operation may be formulated as: 

y(t) =A sin <f> . . (5) 
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Equations (1) to (5) summarize the steps of a tool according to one of the 
described embodiments. When put together, the equations are: 
A = ju x e sin <f>, 
co = ju 2 eA cos <f> , 

5 <j> i = co + /J 3 a>, 

y(t) = A sin <f> , 
e(t) = u(t)-y(t): 

It has been observed that the nonlinear non-autonomous dynamical system 
represented by the above set of differential equations possesses a unique asymptotically 
10 . stable periodic orbit which lies in a neighborhood of the orbit associated with the desired 
component of the function ii(t). In terms of the engineering performance of the tool, this 
indicates that the output of the system y(t) = A sin <f> will approach a sinusoidal component of 
the input signal u(t). Moreover, the slow variations of the parameters in u(t) may be tolerated 
by the tool. v: 

15 The convergence speed in tracking variations in the afnplitude of the input 

signal is observed to be mainly controllable by the assignment of the value of the parameter 
: fi x . The larger p x is chosen, the faster the too! follows variations in amplitude. The inherent 
trade-off existing in tracking capability is with accuracy. The larger ju x is chosen, the higher 
level of error is introduced in the convergence mechanism. Likewise, the assignment of 

20 values of parameters /u 2 and /j 3 provides for controlling phase-frequency tracking speed 

versus accuracy. Roughly speaking, the larger the value of ju 2 is assigned, the faster the tool 
follows variations in the phase and a higher convergence error occurs. 

One way of improving the aforementioned speed/error trade-off is the 
introduction of low pass filters before or after, or equivalently within, the integration 
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operations 20 and 60. This allows for the assignment of large values of parameters /u x and jj 2 
while maintaining the error within a given desired range. 

Figure 1, while outlining the detailed steps of a tool according to one of the 
described embodiments, shows the implementation in the form of the composition of simple 
5 blocks, and therefore presents a tool for the extraction of sinusoids of time-varying nature. 
Such tool may be implemented using digital hardware in which case digital circuits are used 
to perform the required arithmetic operations. An example of such an embodiment of a tool 
may be implemented on a field programmable gate array (FPGA) platform. 

Alternatively, analog circuitry may be employed to construct the blocks of the 
10 tool shown in Figure 1 . • ; 

,■ Other arrangements, such as one employing mechanical components to 
perform arithmetic tasks, may be constructed.. . - 

Yet alternatively, the tool can be easily implemented numerically within: a . > : u? 
software environment. Such a software code may then be executed by the use of computers, 
1 5 microprocessors, microcontrollers or digital signal processor (DSP) platforms or other : % 
computational devices. Numerically, one, and not the only, possible way of expressing the set 
of equations governing the method of the invention in discrete form is: 

A[n + l]= A[n]+ {i x T s e[n]sm0[n\ (6) 
co\n + 1] = co[n] + //^c^cos </>\n\ (7) 
20 ^[« + l]=^[n]+r s <y[n]+// 2 // 3 r ;t e[/i]^[«]cos^[n] > (8) 

y[n]= A[n]swt[n,] , (9) 
e[n] = «[«]- y[n\i (10) 
where a first order approximation for the derivatives is assumed, T s is the sampling time and 
n is the index of iteration. 
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The tool for signal analysis and synthesis is general in the sense that it offers a 
signal processing tool capable of extracting a desired sinusoid which may undergo variations 
in all three parameters: amplitude, phase and frequency. Very often, the component of 
interest is specified in terms of its frequency implying a priori that the frequency is more or 
5 less known and almost fixed. In such situations, the tool implemented on that basis can be 
further simplified. For this matter, supposing that the frequency is more or less fixed around . 
<*>o =, 27 &o 9 one can rewrite the equations for the tool and obtain ' 

•\- A = //,esin^, . 

10 ^ = 0) o + ju 2 eAcos(f>, 

y(t) — A sin \ <f> , 
e(t) = u(tj-y(t). 

It is noteworthy that here the parameter m is reused and is equal to U2U3 35 of Figure 
1 . Another special case is when the parameter wj is taken to be zero, yielding the 
-15 following set of equations: .» 

A = //jesin^, 

&> = ju 2 eA cos 0, 

y(t) = A sin (j> , 

20 e(t) = u(0-y(t). 

The tool outlined by these equations is observed to perform well under more 
or less fixed frequency condition. This is similar to the first special case in terms of 
performance. It is not as simple. The compensating advantage is that like the general case a 
direct estimation of the frequency is available, an (estimation which is more or less a constant 

25 quantity. 
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The two special cases just mentioned do not have to be classified separately as 
they automatically result from the general case under specific conditions. Another special 
case arises when one eliminates the amplitude term in equation (3) of the set of equations. 
This special case is of particular importance in the sense that it does not automatically result 
5 from the general case and has to be considered separately. In this case, the equatiqns may be 
written as m -.■ v. 'f 

• ; - A = /2 x es\x\(f> 9 

d> = ju 2 ecos<fi, 
<fi = co + /i 3 cb 9 

10 y(t) =Xsin <f>, 

: '\ e(t)=u(t)-y(t). 
This alternative procedure and the tool implemented thereupon are also 
observed to perform well under all possible conditions. _ - t V- 

The tool, inathematically described by the equations throughout this 
15 disclosure, refers to a dynamical system which may be described in many other ways. For >r 
example, it is known in the art that a given set of equations can be converted to an alternative 
set of equations by an act of change of variables. For instance, equations (1) to (5), or those 
describing special cases of the latter, are framed in cylindrical coordinate system having 
dimensions {A, <p,a>). Obviously, a change of variables may be employed to frame the same 
20 equations in Cartesian coordinates system. As another example, it is known in the art that the 
parameters //, , /i 2 and , being arbitrary positive niimbers yet to be determined according to 

a specific application, may be replaced by other positive quantities such as // 1 v4 2k , // 2 A 2k and 
M 3 A 2k , A and k being the amplitude and an arbitrary integer, respectively. Variations of this 
sort, exemplified by the two change of variables just mentioned, although result in equations 
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of different appearance, are essentially the same as the general case described in this 
disclosure and are not considered separately. 

The dynamics of the tool presents a notch filter in the sense that it extracts (i.e. 
lets pass through) one specific sinusoidal component and rejects all other components 
including noise. It is adaptive in the sense that the notch filter accommodates variations of the 
characteristics of the desired output over time. c 

The analogy with notch filters should not limit the applicability of the tool. 
From a different perspective a number of core units configured to be employed in parallel, 
can be thought of as a frequency domain analysis tool, such as Fourier transform. From an 
entirely different perspective, the present technique presents a new PLL structure, as will be 
described below. 

To illustrate the performance of the tool, a number of simulations are now ' 
presented. Figure 2 shows the extracted signal of an input of unit amplitude pure. sinusoid / a 
with frequency /- 60 Hz and random phase. The initial conditions are chosen as f 0 = 50 Hz, 

A 0 = 0, and (f> . = 0. With a moderate set of parameters, .//, = 200, ju 2 = 20000, // 3 = 0.02,^ 

the convergence is achieved ; in a few cycles. "» >. 

* • . Figure 3 shows the performance of the topi in tracking time variations in the 
amplitude. A step change of 10% in the amplitude of the input signal, which is otherwise 
taken to be the same as before, occurs. The new value of the amplitude is tracked in a few 
cycles while the phase and frequency of the signal undergo transients for a few cycles. 

Figure 4 shows a similar phenomenon in which the frequency of the input 
signal undergoes a step change. It is observed that the variations are effectively tracked with a 
transient of a few cycles. Values of the parameters are retained the same as before. 

As far as the robustness of the tool with regard to its internal structure, and 
most importantly with regard to the adjustment of parameters, it has been observed that 
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performance of the tool is almost unaffected by parameter variations of as large as 50%. The 
tool is found extraordinarily insensitive to variations of even orders of magnitude in its 
internal setting. This merit is well appreciated in the context of some methods in which 
variations of the order of 0.1% may pose serious stability or performance problems. 
5 In noisy environments, one cannot achieve the maximum speed of 

convergence without a trade-off with accuracy. Figure 5 shows the performance of the tool 
in noisy environments. It is the same scenario as that of Figure 2 with the only difference that 
a white noise at the 20 dB below the level of sinusoidal input (i.e. of 10% of the magnitude) 
is added to the input signal. The extracted frequency is picked as an index of performance. It 
10 is found that the presence of 1 0% noise in the input generates an error of about 0.5% in the 
frequency estimation at the same speed of convergence. It is notable that for each given 
application, one can modify the values of p -parameters to accordingly balance the speed and 

accuracy. •/ 

. The tool is extremely simple. To elaborate further on this fact;, a comparison is 
1 5 now made between the tool and the discrete Fourier transform (DFT). DFT is well-known for ; 
its simplicity of structure and hence its short execution time or little computational demand^if 
implemented numerically in a software environment. Using DFT,. if written out in recursive 
form, the real and imaginary parts of the fundamental phasor, for example, of the input signal 
u(t) can be iteratively computed as: 

20 a\ n ] = a\n-\] + ^{u[n}-u[n-N\cos^), (11) 

b\n] = b\ n -\h^MMn-N]jsm{^ (12) 

where ufn-NJ is the input sample corresponding to the previous period which is saved in 
memory. The fundamental component of u(t) is expressed as: 

u,(t) = A,sm{o) 0 t + e,), (13) 
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in which A, and O x can be calculated from (1 1) and (12) by: 

A x =^lb\, ' (14) 

0 X =tan-'^J. , (15) 

The two recursive equations have very simple structure, Two multiplications 
are the major computations involved. However, as many samples as the window length* 
should be saved in the memory. Note that a relatively large amount of computation is 
required to provide a signal /j x (t) out of a\ and b x . A square root and an arctangent 
operation are to be performed. These two functions are riot straightforward to implement and 
need a considerable amount of memory and hardware. Moreover, further work has to be 
done to synthesize a synchronous signal fu x (t), using a PLL for example. 

Since the tool is not window-based, no input data is required to be saved in the 
memory. It only needs to have the instantaneous value of the input. To compare the tool with 
DFT from the computational order point of view, one needs to consider all the required 
computations by both methods to provide the same result. DFT needs to do all the X 
computations in equations (11) to (15) to provide a fundamental, its amplitude, and its phase 
angle. The tool needs to perform the computations in the set of equations (6) to (10) to J 
provide the same results. It is obvious that the tool lends itself more easily than DFT to the 
practical implementation in terms of computational volume of the software. From a hardware 
point of view, since the extracted fundamental and distortion are directly available, there is no 
need for any synchronizing function such as PLL to implement the tool. 

From one perspective, the tool is an adaptive amplitude, phase and frequency 
estimator 600. This functional description of the invention is pictured in Figure 6, wherein a 
core unit 602 is adapted to receive an input signal 604 and to provide amplitude 606, phase 
608 and frequency 610 output signals. In this sense, the topi belongs to the general category 
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of signal estimation systems and methods such as those widely employed in phasor 
measurement. It can also be thought of as an adaptive notch filter whose center frequency is 
adjustable and is flexibly moving. In this sense, it is comparable with adaptive filters. It can 
be envisaged as a signal analysis topi which analyzes individual components of frequency 
5 decomposition of a given signal. In this sense, it is comparable with DFT. 

Unlike most signal processing techniques which provide estimates of 
characteristics , of the components of a given signal, the tool not only does provide estimates 
of the amplitude, total or constant phase and frequency, but also synthesizes the desired 
sinusoidal component itself in real time. This is feasible due to the fact that the total phase is 
10 estimated at each instant. From this perspective, the present invention is a new PLL 700 as 
pictured in Figure 7, wherein a core unit 702 is coupled to receive an input 704 and provide a 

'>*.>' K 

synchronous output 706. , . 

The input signal, out of which the system is supposed to extract a clean single 
sinusoid, may contain components lying in frequency bands far off that of interest. A filter 

1 5 may be inserted at the input of the core unit to eliminate some of the undesired components 
and thus enhance the speed-accuracy tradeoff of the tool. Obviously, this filter does riot need 
to be sharp meaning that its structure is very simple and the phase delay will not be s * >*■ 
excessively long. The phase delay and the gain of the filter will be functions of frequency. 
The tool, by the very fact that it synchronously estimates frequency, amplitude and total 

20 phase, allows for the correction of the phase delay and change in amplitude. The estimation 
of frequency at each instant may be used to obtain the phase delay caused by the filter. Such a 
quantity is then subtracted from the total phase which is also instantly generated. Likewise, 
the amplitude can be corrected. The desired output can then be simply synthesized as 
illustrated in Figure 8. As shown in Figure 8, an input signal 804 is coupled via a filter 806 to 

25 a core unit 802. The core unit 802 provides amplitude 808, phase 810 and frequency 812 
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outputs. The outputs 808-812 are coupled to an amplitude and phase correction unit 814 
which provides an amplitude and phase corrected output 816. The sine of the corrected 
output is provided by a sine component 818, and the product of the sine and the output 816 is 
provided by the multiplication function 820 to provide an output signal 822. A particular 

v 5 example of an application of the configuration of Figure 8 is a power signal interference 
eliminator. Such a device may include a filter such as that shown in Figure 8 to extract the 
interfering po>yer line signal that is then subtracted from the input to provide a clean signal 
out of which the power signal is removed. 

As a signal analysis tool, a number of core units, or configurations of Figure 8 
, .10: as core units for improved performance, may be connected in parallel, or cascaded, to 

decompose the multi-component input signal into its constituting sinusoidal components. 
Figures 9a and 9b shows two possible combinations 900 and 900\ In Figure 9a, a plurality 
of core units 902 are coupled in parallel to receive an input signal 904 and to respectively ^ 
provide output signals 906-910. In Figure 9b, a plurality of core units 912-916 are couplecTin 

15 series/cascade. . The first core unit 912 is coupled to an input signal 918 and provides an 

output signal 920 and error signal 922. The error signal 922 is coupled as an input to the/core r 
unit 914, which provides an output signal 924 and an error signal 926. The error signal 926 is 
coupled as an input to the core unit 916, which provides an output signal 928. 

Of course, one can conceive of other combinations to build multi-core units 

20 tailored by needs of any given application. For example, if the desired component is masked 
by components of a more or less known frequency composition, undesirable components may 
be extracted and subtracted from the input signal first to provide a cleaner input for the core 
unit which is supposed to extract the desired component. Figure 10 shows such a 
configuration. In Figure 10, a plurality of core units 1002-1006 are coupled to received an 

25 input signal 1008 and to provide corresponding output signals 1010-1014. The output signals 
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1010-1014 are subtracted from the input signal 1008 by a subtractor 1016, the output 1018 of 
. which is coupled as an input to the core unit T020. The core unit 1020 provides an output 
1022 based upon the signal 1020. It is noteworthy that each core unit tries to extract a ; 
sinusoidal component of the frequency that is closest to its specified initial condition. In order 
5 to avoid overlapping duties of the core units, the operating frequency of each of the core units 
has to be within a specified range. This can be achieved by introducing range limiters within 
frequency integration 60 of Figure 1, for instance. 

It will be appreciated from the foregoing discussion that the tool herein 
described, which is well suitedTor signal extraction and particularly of low level signals in 

10 noisy environments, may be adapted for use in audiometry to improve sensitivity and to. Y 
reduce testing time. The following discussion:and examples therefore focus on application of 
the tool in audiometry. \c 
Referring to Figure 11, an audiometer 1 100 includes three main modules: the- *g 
data acquisition/transducers module 1 1 02, the signal processing module 1 1 04 and the user ' 

1 5, interface 1 1 06. A probe 1 1 08 is coupled to the data acquisition module 1 1 02, and a personal # 
computer or other computational/data management device 1 1 10 may be coupled to the signal 
processing unit: The acquisition/transducer module 1 102 and the signal processing module 
1 104 are described in more detail below. The user interface 1 106 may include a data entry 
device, such as a keypad, keyboard, touch screen, mouse or other pointer device, voice 

20 activation and the like for accepting commands, data input and actions from a user. The user 
interface 1 1 06 may also include a display and/or printer device for providing information to 
the user. The particular configuration of the user interface is not critical to operation of the 
audiometer 1 100, and the user interface function may also be provided by the personal 
computer 1110. Thus, one of skill in the art given the nature of the data and commands to be 

25 accept and the information to be provided to the user would be capable of specifying a 
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suitable user interface. For example, the user interface may be a laptop or palm-type 
computing device with a graphic iiser interface and a keypad, mouse, joy-stick, touch pad 
and/or voice recognition input device. 

Data acquisition unit 1 1 02 is the medium between the processing unit 1 1 04 
and the probe 1 108 which transmits and receives signals to the testing subject. For example, 
in DPO AE; the probe 1 108 may include a transducer 1202 and a microphone 1204 for 
transmitting and receiving signals in the human audible range as shown in Figure 12. For 
SSAEP, auditory brainstem response (ABR) or other types of evoked response audiomery, 
the probe may include an audio transducer for introducing a stimulus to the testing subject 
and electrodes for detecting an electroencephalogram (EEG) of the testing subject. It should 
be clear, however, that the probe 1 108 may be configured for prbviding virtually any type of 
stimulus: auditory, electrical, haptic and/or visual and for receiving virtually any type of ^ 
response signal: auditory,. electrical; haptic and/or visual in connection with audiometry f ^; 
and/or other forms of objective patient testing. In that regard, one of the main functions of? 
the module 1104 is to convert using a digital-to-analog converter 1206 digital, signals V 
produced by : the signal processing module 1102 to stimlus signals, typically analog signals^ 
which are then -conditioned by optional filtering 1 208 and amplification 1210: and converted 
to audio signals for transduction to the testing subject by the transducer 1202. Conversely, the 
response signals, audio or electrical signals as the case may by, recorded by the probe, are 
conditioned by j .amplification 1212 and filtering 1214 and converted be analog-to-digital 
converter 1216 to digital signals to be processed by the signal processing module 1 102. 

The signal processing module 1 102 of the audiometer 1 100 produces the 
digital form of the artifacts and extracts and measures the response signal. A DSP, or if the 
computational/architectural demand is low even a microcontroller, can be employed as the 
hardware platform of this unit. Signal processing is embedded as the software in such a 
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hardware platform. Alternatively, and provided that the complexity of the signal processing 

; algorithms remains low, signal processing unit may be implemented solely in hardware using 
programmable logic array (PLA) or field programmable gate array (FPGA) technology. In an 
ideal case, namely when the signal processing algorithm is not excessively complex, the 
5 hardware does not require a PC for its operation; however, interfacing to a PC is usually 
provisioned for data management. ' . 

Figure 13 shows the main functions of the software embedded in the signal 
processing module 1 102. The software is essentially responsible for the generation of the 
artifact signals and extraction of response signals, i.e., processing portion 1302, as well as 

10 management of the input/output of data 1304-1308. 

; With reference now to Figure 14, a signal processing scheme for the extraction 
of DPOAE is illustrated. One of ordinary skill in the art will appreciate that the principles;-/ 
herein described for the configuration of an audiometer, and particularly for extracting a \i£ :> 
response signal from artifacts and noise, may be used to construct a signal processing scheme 

15 for extracting transient or steady-state evoked response signals, other forms of otoacoustic&r 
emission signals,and the like. The following example of a single processing scheme adapted 

i for DPOAE therefore illustrates the general principles of the invention herein described arid 
its adaptability and in no way should be considered limiting of the applicability of the general 
signal processing scheme to a particular form or type of audiometry. Moreover, it will be 

20 appreciated that it may be possible and/or desirable to provide input signal filtering, such as 
described in connection with the embodiments of a core unit illustrated in Figure 8, and/or to 
provide enhanced core units, such as those illustrated in Figure 10 ; to enhance signal 
extraction. Also, signal sampling frequency may be modified, i.e., increased or decreased, 
either alone or in combination with input signal conditioning techniques to enhance signal 
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extraction. These input signal conditioning and sampling techniques may be particularly 
advantageous when attempting to extract auditory evoked potential signals. 

The scheme employs three core units 1402-1406 toxonstruct a high 
performance DPOAE extraction module 1400. Each core unit 1402-1406 is capable of 
focusing on and extracting a pre-specified sinusoidal component of its input signal which 
may contain many other components including noise. More importantly, each can effectively 
follow variations in the amplitude, phase (and frequency) of the extracted sinusoidal 
component. The underlying mathematics ensuring stability and performance of such core 
units 1402- 1406 may be in accordance with any of the foregoing described structures, and 
thus,"the structure of the core units remain extremely simple. They are found to be very ,\ 
robust with respect to variations in the internal settings as well as external conditions arid 
exhibit superior performance over. existing linear adaptive and FFT-based algorithms. >..'■>... 

sThe input signal is often assumed to consist of two pure sinusoids with Xry t 
frequencies /, v and f 2 at a very high level (usually about 60 — 70 dB) and a . very low level . 
DPOAE 2 f } -]f 2 at about 0 dB. It is polluted by a noise usually considered to be at about--T0 
dB level. The noise in fact represents the totality of all undesirable signals that may be : \:y 
present in the environment in which the examination is being conducted as well as 
unavbidable white Gaussian noise. While it is possible to use a single core unit to extract the 
DPOAE signal, because of the degree of pollution (artifacts and noise) use of multiple core 
units may provide better performance. 

The first two core units 1402 and 1404 are set to extract the artifacts of the 
stimulus signal, i.e., frequencies /, and f 2 . They effectively do so with very small errors. 
The extracted artifacts are then subtracted from the input to produce a signal 1408. DPOAE 
now has a higher relative portion of the signal 1408 because the artifacts have now been 
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removed. The third core unit 1406 is then set to extract DPOAE. That is, the third core unit 
1 406 is configured to extract the very low level DPOAE 2f v -f 2 signal. 

The core uhits 1402-1406 may be configured in digital or analog hardware or 
software with good results. The core units may also be configured in accordance with any of 

5 the algorithms herein described for the extraction of signals. 

< 

The performance of the audiometer 11 00 configured to extract DPOAE can be 
measured with the use of computer simulations as demonstrated by the following examples, 
which are intended to demonstrate the capability of the audiometer 1 100 and should not be 
takeri as limiting of the scope of the patent. - 

10 The arrangement of Figure 14 including the three core units 1402-1406 is used 

to extract DPOAE in a : simulated environment All the initial conditions are taken to be zero! 
The first two units 1402 "and 1404 are employed to extract the artifacts and must do so with as 
small an error as possible to produce an input to the third unit containing DPOAE signal as. a 

* main component. Otherwise, the third unit may not be able to extract DPOAE signal with v 

15 sufficient accuracy. Therefore, the parameters may be chosen relatively small, and for tHe| * 
described example they riiay be chosen as: p x = 100 and ju 2 = 50. With these values, it is ; 
expected that the convergence is achieved within about 1 00 ms. The same setting may also be 
used for the third unit. Two forces oppose each other in the selection criteria for the values of 
the parameters of the third unit. On the one hand, since SNR is very low for the third unit 

20 (DPOAE being taken as the signal and the total residual errors of the first two units as the 
noise), the parameters must be very small. On the other hand, since the input has a very low 
level, the value of /u 2 must be large. These two forces tend to cancel each other for the case of 
ju 2 so that the original value of 50 remains a good choice. As for the value of ju x , one 
expects that it should be small to achieve acceptable accuracy in the extracted signal. This is 

25 true, but since only the amplitude of the extracted signal (i.e. the level of DPOAE) is of 
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interest, one can substantially increase //, . The accuracy in the estimated extracted signal 
(and not necessarily its amplitude) is thus traded off, which in this case is not of importance. 
The estimation of the amplitude remains accurate while the convergence speed is high 
enough. A rough setting of the values of parameters /i, and // 2 can be done. The core unit 
algorithms have been shown to be very insensitive with regard to the settings of the 
parameters; variations as much as 50% in the values of these parameters seem to have 
practically no effect on the performance. The Value of the amplitude directly estimated by the 
third unit was used as the estimate of the level of DPOAE. 

The input signal is synthesized by the summation of two artifacts at 
frequencies /j ' = 3000 Hz and f^= \^2 f x = 3600 Hz and the DPOAE signal at 2 /; - / 2 = 

2400 Hz. The amplitude of the artifacts are taken as 1 v and the DPOAE level is 0.5 mv: 

■'*.' t- ' • . - : , % • ■ . • ... 

With these values, the level of the artifacts is set at about 66 dB higher than that of DPOAE. 

Figure 15 shows the performance of the audiometer 1100. It is confirmed that the 

-,• ' ■ " . . • • '.■ 

measurement is completed in about 1 00 ms. 

To present the dynamic performance of the audiometer 1 100, the response of " 
the algorithm to a step change in the level of DPOAE signal is presented in Figure 16. It is' , 
clear that the algorithm follows the Variations in the level of DPOAE at a sufficiently high 
speed. : ' 

A number of measurements can be made consecutively. Figure 17 shows such 
an example. In each case, the frequency of the extracted DPOAE signal is shown. As the 
frequency decreases, less number of cycles are available in a given time interval. For 
example, there are 160 cycles available in 100 ms of a signal of frequency / = 1600 Hz 
whereas in the case of a signal of frequency / = 4000 Hz the available number of cycles are 
400 in the same time interval. This explains why the estimation of DPOAE level is better 
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achieved at higher frequencies. It should be noted that within the range of frequencies of 
interest for DPOAEs, the algorithm provides acceptable estimates in about 100 — 200 ms. 

As noted before, one of the main features of the employed core algorithms is 
their noise immunity. Therefore, one expects that the overall DPOAE estimation scheme be 
immune to background noise. Usually, otoacoustic tests are conducted in quiet examination 
rooms so that the level of noise floor is low enough for good results. Under such controlled 
conditions, a noise floor of about 10 dB lower than DPOAE signal (SNR=10 dB) exists. It is 
obvious that an algorithm which is able to extract a DPOAE signal.highly immersed in 
background noise is desirable in that it alleviates the need for special quiet examination 
facilities such as sound-proof booths. Figures 18a and 18b show the performance of the 
audiometer 1 100 in a highly noisy environment where the level of the DPOAE signal is set to 
be at the level of the noise floor (SNR=0 dB). The incurred error is negligible. In Figures 19a 
and 19b, the level of noise is elevated to be 10 times larger than that of DPOAE signabitself. 
The incurred estimation error is less than ±15% which for most practical application is ; 
sufficiently good. The excellent noise immunity of the audiometer 1 100 not only obviates the 
need for sound-proof examination rooms, but also provides the way to reduce the level of the 
artifacts for more patient- friendly tests. 

The invention has been described in terms of several preferred embodiments 
and examples. One of skill in the art will appreciate that the invention may be otherwise 
embodied without departing from its fair scope, which is set forth in the subjoined claims. 
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